function out = mapFeature(X1, X2)
    %   Degree of the mapping
    degree = 6;
    %   Initialize the result
    out = ones(size(X1, 1), 1);
    for i = 1:degree
        for j = 0:i
            %   Mapping the inputed data
            out(:, end + 1) = (X1.^(i - j)) .* (X2.^j);
        end
    end
end
